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Abstract 

Constants of motion, Lagrangians and Hamiltonians admitted by a family of relevant nonlinear 
oscillators are derived using a geometric formalism. The theory of the Jacobi last multiplier allows 
us to find Lagrangian descriptions and constants of the motion. An application of the jet bundle 
formulation of symmetries of differential equations is presented in the second part of the paper. After 
a short review of the general formalism, the particular case of non-local symmetries is studied in 
detail by making use of an extended formalism. The theory is related to some results previously 
obtained by Krasil’shchi, Vinogradov and coworkers. Finally the existence of non-local symmetries 
for such two nonlinear oscillators is proved. 
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1 Introduction 

The interest of differential geometric techniques in the analysis of systems of ordinary differential equa¬ 
tions (ODEs) is undeniable. Lie point symmetries, integrating factors and their generalizations are just 
several examples of geometric tools which have been successfully applied in the study of systems of ODEs 
and their related mathematical and physical problems. 

In this work, we survey the geometric theory of Jacobi multipliers mm and non-local symmetries mm 
to study a family of relevant nonlinear oscillators that have been attracting some attention in recent years 
[3 nil [Hill [To] [n]. For instance, it was proved that they can be understood as oscillators in manifolds of 
constant curvature m, they admit compatible bi-Hamiltonian structures, and their properties can also 
be analyzed by means of coalgebra techniques [13 E]. Some of their properties have also been obtained 
by means of the so-called A-symmetries [S]- 

First, we use Jacobi multipliers mm to go over the above-mentioned oscillators from a geometrical 
point of view. This allows us to obtain some of their constants of motion. We also obtain known and 
new Lagrangian and Hamiltonian functions for these oscillators. It is worth noting that the new derived 
Lagrangian functions are of a non-mechanical type, i.e. they do not possess a kinetic term given by a 
2-contravariant tensor field. 

Subsequently, we review a ‘new method’ to obtain non-local symmetries developed by Gandarias 
and coworkers 13 [g [7]. We show that their procedure is a consequence of the non-local symmetry idea 
formalised by Krasil’shchik and Vinogradov several years before [I3II1II3. Despite this, Gandarias 
and coworkers’ applications of this method are still relevant, as they illustrate that certain systems of 
differential equations with no classical point symmetries can admit non-local symmetries that lead to 
unveil their properties. 

The study of non-local symmetries demands the use of infinite-dimensional jet bundles [16]. We il¬ 
lustrate how this geometric approach can easily be applied to study our family of oscillators. Indeed, 
the calculation of such non-local symmetries in the problem under study is very similar to the case of a 
finite-dimensional jet bundle J^tt. It is essentially the geometrical interpretation what differs. Addition¬ 
ally, our techniques provide very simple and relevant examples of finite-dimensional diffieties describing 
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ODEs, which is not the usual approach as they are mainly concerned with infinite-dimensional manifolds 
describing PDEs. 

The use of infinite dimensional jet bundles provides other advantages. Many structures of J^tt, e.g. the 
Cartan distribution, become simpler when defined on the infinite-dimensional jet bundle Moreover, 

we can define geometric structures on this latter space that cannot properly be defined on J^tt, e.g. the 
total derivative. Moreover, J°°Tr possesses a geometric structure richer than that of J^tt, e.g. it admits 
Lie symmetries than cannot be described in terms of Lie symmetries defined on p-order jet bundles m- 

Apart from showing the usefulness of infinite-dimensional jet bundles and reviewing Gandarias’ re¬ 
sults, we devise a new idea to easily determine non-local symmetries for certain systems of higher-order 
differential equations. As an application, we retrieve in a natural and rigorous way a result given by 
Gandarias as an ansatz for the oscillators of this work . 

The paper goes as follows. Section 2 is devoted to surveying the theory of Jacobi multipliers and 
its relation with Lagrangians and constants of motion. In Section 3 we apply Jacobi multipliers to 
study relevant types of oscillators. In Section 4, non-local symmetries of differential equations are briefly 
presented. We show that the method developed by Gandarias and coworkers 016 ] reduces to the non-local 
symmetry concept developed by Krasil’hinski and Vinogradov and we apply this idea to the mentioned 
oscillators in Section 5. A method to improve the derivation of such non-local symmetries is provided in 
Section 6 and we summarise our results and comment on our future work in Section 7. 


2 Jacobi multipliers, Lagrangians and constants of motion 

Let M stand for an oriented manifold, i.e. M is equipped with a volume form fl. Given a vector field X 
on M, we call divergence of X relative to ft the unique function div V : M —> R satisfying mm- 

Cxft = (divAT) ft. 

A Jacobi multiplier for Af is a non-vanishing function p : M —>■ R satisfying that p{i{X)ft) is a closed 
form, or equivalently CxitJ-ft) = 0. In other words, p, is such that 

div(/iAr) = 0, p{x) 7 ^ 0, Wx G M. 

From Cxifft) = {Xf + f div X) 12, for every / G we see that a function p is a Jacobi multiplier 

for AT if and only if p does not vanish and satisfies 


pdivX -G Xp = 0 . 


( 2 . 1 ) 


A function / : M —>■ R is a first-integral of X, i.e. XI = 0, if and only if dl annihilates the generalized 
distribution V = {v G TM \ dp G M, v = Xp} generated by X, i.e. (d/)p(Xp) = 0 at each p G M. Let us 
restrict ourselves to an open subset U = {p G M Xp ^ 0 , (d/)p 7 ^ 0} of a two-dimensional manifold M. 
In this case, i{X)ft defines a one-dimensional codistribution annihilating U and hence X admits a Jacobi 
multiplier p such that pi{X)ft = dl. Using that (f(X)U) A d/ = —{i{X)df) Aft = —{Xf)ft, for each 
/ G C°°{U), we find that dl A df = —p {Xf) ft. This expression also shows that / is a first-integral of 
X if and only if dl A df = 0, i.e. / is a function of / and / = (p{I) for a certain real function (p : R —R. 

Jacobi multipliers have many applications [HlIIlIIllIISlEolEIlEg. In particular, we are interested 
in their use to construct Lagrangians and constants of motion for second-order differential equations 
[niMiiiniiis]- Let us briefly survey this topic. 

Assume hereafter that TR is endowed with the volume form U = dx A du. Consider a second-order 
differential equation 


d^a; ^ / dxN 

lW = ^G-di)' 


X G 


( 2 . 2 ) 
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with F : TR —?► R being an arbitrary function. Adding a new variable v = dx/dt, we see that (12.21) can 
be rewritten as 

dec 

TIT = 

f (2.3) 

whose particular solutions are integral curves of the vector field on TR given by 


d d 

r = ,_ + F(*,„)_ 


(2.4) 


A first-integral of T is usually called a constant of motion for T or, equivalently, for system (ESI). By 
substituting v by dx/dt, this first-integral gives rise to a constant of motion to (12.21) . As divT = dF/dv 
in this case, then the Jacobi’s multiplier condition (EU amounts to 


ox ov 


(2.5) 


The Jacobi multiplier ^ satisfying this condition is also called a Jacobi multiplier for the second-order 
differential equation (12.21) . 

Theorem 2.1. The differential equation determining the solutions of the Euler-Lagrange equation defined 
by a regular Lagrangian function L(x,v) possesses, when written in its normal form a Jacobi 

multiplier given by the function 

d^L 


T = 


dv^' 


( 2 . 6 ) 


Conversely, if p is a Jacobi multiplier for a second-order differential equation h2.SM . then admits a 

regular Lagrangian L{x,v) satisfying i t A dl) . 

Proof. Assume L to be a regular Lagrangian for (12.21) and define the non-vanishing function p by (12.61) . 
Note that then F is given by 

1 fdL d^L \ 
p \dx 




and, using this, we see that 


dp d 


+ — (pF) = V— -\- 


dp d fdL d^L \ d^L d^L 


dx dv 


dx dv V dx 


— V 


dxdv J 


= V 


d^L d^L 

— u „ „ = 0. 


dv^dx dxdv dxdv dv'^dx 


Therefore, p given by (12.61) satisfies the Jacobi multiplier equation (12.51) . Since L is assumed to be regular, 
then the function p given by (12.61) does not vanish and becomes a Jacobi multiplier of (12.21) . 

Conversely, if /r is a Jacobi multiplier for (12.41) . then the functions L satisfying (|2.6I) are of the form 


/ V nv' 

dv' p{xX)d((j)i{x)v(j>2{x) 


(2.7) 


for arbitrary functions ((ii, <('2 : R —>■ R. The term (pi (x) u is a gauge term which can be hxed equal to 
zero, i.e. L can be assumed to be of the form 


/ V nv' 

dv' p{x,()dC-\-(p2{x). 


( 2 . 8 ) 


and the function <p 2 can be chosen in a unique wav|27J. up to a constant, so that the Euler-Lagrange 
equation reproduces the equation for the integral curves for the given vector field p.4|) . Indeed, using 
we see that 


dL r , r dp ,,,,^dp2, , 


3 









and in order to the Euler-Lagrange equation for the Lagrangian (I2.8|) to give the dynamics we should 


have: 






d Li 

-^{x, C) dC + fj.{x, v)F{x, v). 


But note that 


dv 



dx 


ix,0 dC + ^i{x,v)F{x,v) 



dx 


(a;,C)dC 


dfi dF 

= + Pkx,v)^{x,v) + ^i[x,v) — {x,v) 


(2.9) 


which vanishes because of the multiplier condition (12.51) . Consequently, the function (j )2 exists and is 
uniquely determined, up to a constant, by 


d(j)2 

dx 


(x) = ' 


d^jL 

dx 


(x,C) dC + ^J.{x,v)F{x,v) 



dfi 

dx 


(x,C) dC- 


( 2 . 10 ) 


An integration by parts in the double integral leads to 




d^ 

dx 


(x,C) dC = 



9/r 

dx 


(x,C) dC 


v' |^(x,u') dx'. 


that when substituted in (12.101) gives 

^(x) =/x(x,x)F(x,x)+y c|^(a;,C)dC, 
which using the Jacobi multiplier condition (j2.5l) reduces to 


d(j)2 

dx 


(x) = fj.{x, x)F(x, v)- f {x, C) dC, 

Jvn oC 


and therefore (j )2 is given by: 


(l)2ix) = J (/rE)(C,xo)dC. (2.11) 

Additionally, since /i is non-vanishing and in view of (12.6|) . we obtain that L is regular. □ 

Another remarkable result concerning the inverse problem is given in the following Proposition [281129] . 

Proposition 2.2. If I is a constant of motion for the vector field T at a point ^ G TR where Pj 0 and 
(d/)^ 0, then 

L{x,v) = vJ (2.12) 

is a Lagrangian for the given vector field around a neighborhood of 


Proof. Since P and dl do not vanish at there exists around this point a Jacobi multiplier p for P 
relative to fl = dx A dx such that 


pi{r)Ll = dI <;=^ p(v dv — F{x,v) dx) = dl. 


Therefore, 


1 ^ 
X dv' 


-pF = 


dx 


4 






In view of Theorem 12.11 there exists a Lagrangian L such that 


dv"^ 


1 dl 
V dv’ 


from where, an integration by parts leads to 
dL _ r 1 dl _ I{x,v) 


/ 


d I 


H^,0 


dc 


This shows that L must be given by 


L{x,v) = vJ '^^^ 2 ^^ dC + a(x) 


for a certain cj) : x (j){x) G K. Imposing L to be a Lagrangian for T and recalling that dl/dx = —fiF, 

we obtain 


dL d dL ^ 
dx~ dtdv ~ ^ 


f 


[liF]{xX) ^ f 

da; dt \ V 


f 


Hx,0 


dC = 0. 


Hence, 




Fdl 

V dv 


F 


[^^F]{x, c) _EL = ^ = Q 


e 


dx 


Hence, (j) is an irrelevant constant and we obtain that L is essentially given, up to an also irrelevant gauge 
term, by (12.121) . □ 

Apart from providing a variational description for second-order differential equations (SODEs) as 
(1^ . Jacobi multipliers can also be employed to derive their t-independent constants of motion, namely 
first-integrals for the associated T. More specifically, given two Jacobi multipliers fii and fi 2 of the vector 
field r, the function 

Ml 

= — 

M2 

is a constant of motion for (12.31) and, by substituting v by dx/dt, we obtain a t-independent constant of 
motion for (12.21) . Indeed, as div(^ir) = T^i + fii divT = 0, for z = 1, 2, it follows 

^ ^ fl2^fll - flirfi2 ^ g 

mI 

Consequently, the non-uniqueness of such a Lagrangian function, i.e. the existence of alternative 
Lagrangians can be used to determine constants of the motion as it was proved in |30j for the one¬ 
dimensional case and generalized in [3T] for the multidimensional case (see also [32] for a geometric 
approach). 

In addition, given a non-vanishing t-independent constant of motion ip for T, then piip is a new Jacobi 
multiplier for T. This shows that given a fixed Jacobi multiplier pi for T, any other Jacobi multiplier fi 
for r arises as the product of /ii times a non-vanishing function G G C°° (K.) of a given nontrivial constant 
of motion ipi for T, i.e. fi = G{ipi)fii. 

3 Jacobi multipliers and nonlinear oscillators 

Let us now use the above results to analyse the nonlinear oscillators 


d X 

dF 


kx 


1 -I- kx'^ 


dx Y 
dt) 


a^x 
1 -I- kx'^ 


= 0 , 


a G 


(3.1) 
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and 


da; 

dt 


2 


a G K, 


(3.2) 


d^a: 

df^ 


kx / 
1 + kx'^ \ 


+ 


a'^x 


(1 + kx^y 


= 0 , 


which have recently been drawing some attention [3 ESI [El Ea [To] . For instance, the Hamiltonian 
description of the quantum analogues of both systems led to suggest a Lagrangian description for nonlinear 
oscillators [33] . Here, a; G K when k > 0 but |a;| 7 ^ 1/^/\k\ for k < 0. For simplicity, we study the bounded 
motions with x G (—l/\/—k, l/\/—k) when fc < 0. Some generalizations of these results to higher-order 
dimensions for (EH) were devised in [12] and some of the non-local symmetries for EH and EH were 
described in [5]- 

The second oscillator EH is the one-dimensional case of the Hamiltonian superintegrable system 
studied in m- Apart from its superintegrability, this system has attracted some attention due to the 
fact that it can be investigated through an s[(2, ]R)-Poisson coalgebra (cf. [IH])- The straightforward 
generalization to n-dimensions of both oscillators leads to oscillators of constant and variable curvature 

m- 


3.1 First nonlinear oscillator 

We can write EH as a first-order system: 

{ dx 
dt 

du — a‘^)x 

dt 1 + kx'^ 

System (13.31) describes the integral curves of the vector field 

d x{kv'^ — a^) d 

^ dx 1 -I- kx"^ dv ’ 


(3.3) 


and as divF = 2kxv/{l + /cx^), its Jacobi multipliers, fx : TR —>• M, are the non-vanishing solutions of the 
equation 

div {fjT) = F/i -|- /idivF = 0, 

that in this case can be written as 

dfj. x{kv‘^ - a'^) dfi 2^ikxv 

+ i + fa^ ^ 

This equation can explicitly be solved by the method of characteristics. This method reduces solving the 
above PDF to determining the integral curves of a vector field: the so-called characteristic curves. When 
characteristics are determined, solutions for the PDE are obtained by gluing them together giving rise to 
a hypersurface (see [34] [35] for details). The characteristic curves of p.4p can be described by means of 
the referred to as characteristic system [341135] , namely 


dx (1 -I- kx'^)dv (1 -I- fcx^)d/r 
V x{kv'^ — a^) 2kxvyL 


Let us solve these equations for k = 0, i.e. the harmonic oscillator. In this case, we have divF = 0 and 
Jacobi multipliers become mere first-integrals of F. The characteristic curves are given by 


/r = Ti, x'^ = T2, 


where Ti,T 2 are real constants. This means that Jacobi multipliers are non-vanishing functions of the 
form /r = /r(x^ -|- u^). 
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By integrating the characteristic equations for fc ^ 0, we find that the characteristic curves are given 


by 

where Ti, T 2 are real constants. We know that any surface (a;, y, y,{x, v)) obtained by gluing characteristic 
curves is a solution to (13.41) . namely any subset (x, y, y) of points of satisfying 

K{{1 + kx^)/{kv'^ — a^), /r(l + kx^)) = 0, 


for a fixed function if : —>■ R, with d^K ^ 0 (observe that this amounts to ^ = /i(x, y)). If we impose 
the boundary conditions ^i(0,u) = 1 and ^ 2 ( 0 ,u) = l/(fcu^ — a^), we easily obtain, respectively, the 
Jacobi multipliers 

_ 1 _ 1 
1 + kx'^ ’ kv"^ — 

Both Jacobi multipliers lead to the existence of a constant of motion for T of the form 


^ _ M 2 _ 1 + kx^ 

Ml kv'^ — 


(3.6) 


Note that, as previously remarked, a constant of motion I for (12.21) and a Jacobi multiplier mi enable us 
to recover all the Jacobi multipliers for X as y, = G{I)yi with G(i) being an appropriately non-vanishing 
function. 

Let us now turn to working out a Lagrangian for oscillator (13.ip by using the method of Jacobi 
multipliers and mi (observe that mi is also a Jacobi multiplier of (13.ip for fc = 0). Recall that this method 
states that dan admits a Lagrangian Li : TM —>■ R satisfying 


d^Li _ _ 1 

dv^ 1 - 1 - kx'^ 


This yields 

1 

Li{x,v) = 2 (i + fca;^) + 4>i{x) v + 4>2{x), 

for certain functions 4>i^4>2 : R —>■ R. We can set the gauge term (j)i{x) v equal to zero while (j )2 is to be 
determined using (12.lip . More specifically, choosing vq = 0 in Theorem 12.11 we obtain that the function 
02 is given, up to the addition of a constant, by 


02 (a:) 


'(MF)(c,o)dc= r 

Jo 


1 (-a 2 )C 


1 -k 1 -k kC 


dC 


r 2 1 ~\^ 

a 1 

_2fc 1 -k kC^ _ g 


o? X 


2 


2(l-kfcx2) ■ 


Therefore the Lagrangian for (13.11) is given, up to addition of a gauge term, by 


Li{x,v) 


Iv^ — a^x^ 

2 1 -k kx'^ ’ 


and the corresponding momentum and Hamiltonian read 


dLi V 

^ dx 1 -k kx'^ 


1 


= 0(1 + kx )p -k - 


1 


2 1 -k fcx2 ■ 


We note that the function Li coincides with the Lagrangian obtained in |33j by direct approach. It is 
remarkable that Li is a standard mechanical Lagrangian: it is given by a kinetic term quadratic in the 
velocities minus a potential term. Moreover, we can also prove that Li is, up to an irrelevant additive 
constant, the Lagrangian function L for T obtained by using Proposition 12.21 and the constant of motion 
/ = Mi/( 2 fcM 2 ): 


L(x, v) = v J —^ 5 —dC = v J 


kC,"^ — o? 


2 fc(l -kfcx 2 )C' 


rdC = 


kv^ 


2 fc(l -k fcx 2 ) 


— Li 


a 
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Meanwhile, the second Jacobi multiplier, ^,2 = gives rise to a non-mechanical Lagrangian. 

Indeed, if we assume k > 0 and d^L 2 ldv^ = /i 2 , we obtain that, up to irrelevant gauge terms, the 
corresponding Lagrangian reads 


L2(x,v) = -^^arcth 

yka 


V , (Vkv\ 1 , , 2 ,2 


I - - kv^\+(l)2ix), 

a 2k 


where 


Hence, 


Mx) = [ (M-F)(C,0)dC = f 

Jo Jo 


CdC ln(l -I- kx'^) 


0 l + fcC^ 


2k 


L2{x,v) = -r^arcth 

Vka 


\fkv \ 1 1 + kx 


in ■ 


In consequence. 


p= -- 


arcth(-\/fcv/a) 
\/k a 

The case fc < 0 leads to a similar result. 


H2{x,p) = 


a j 2k — fcu^| 

I — th(-\/fcpa)^ 


I J- kx"^ 


3.2 Second nonlinear oscillator 

We now apply Jacobi multipliers to nonlinear oscillators (132]). Proceeding as before, we consider such 
systems as first-order systems by adding a new variable u = dx/dt to obtain 


= V, 


kxv^ 


1 J- (I -I- kx'^Y 


(3.7) 


We drop the case fc = 0 as it leads to the standard harmonic oscillator which was analysed in previous 
section. So, we now focus upon the case fc 0. The multiplier equation (12.51) for the vector held associated 
to the above system reads 


/ kxv^ a^x \ dn 2iJ,kxv 

^ dx \1 + kx'^ {1 + kx'^)^ ) dv 1 + kx'^ 

whose characteristic system is 

dx (l-l-fcx^)^du (l-l-fcx^)d^ 

V kxv"^ (\ + kx'^Y + X 2kxvyL 

The equality between the hrst and the last term shows that 

dx ^ d^ ^ fj. ^ 

I -I- fcx^ 2kxyi 1 -I- fcx^ 

for a real constant Ti. Meanwhile, the integration of the equation 

dx (l-l-fcx^)^dv 

kxv^{l + kx'^y- + o?x 


V 
























goes as follows. Rewrite the above equation as 

kv'^ {1 + udu dv"^ 

{1 + kx^)^ xdx dcc^ 

The local change of variables w = , z = x'^ transforms the above equation into 


drc k o? 

dz 1 + kz (1 + kz)^ ’ 


which finally gives 


fc(l + kx^)^v^ — 

k{l + 


for a certain real constant T 2 . Resuming, we obtain the characteristic curves 




1 + kx'^ 


= Ti 


k{l + kx^fv^ — a 

fc(l + kx"^) 


= T 2 


Imposing, for instance, /i(0,u) = 1 or ^j.{0,v) = {kv^ — o?)lk^ we obtain 


^1 = 1 + kx^, fi2 = ^(1 + kx^Yv^ — o? , 


which enable us to define a constant of motion 


/r2 fc(l + kx^)'^v^ — o? 

fii 1 + kx"^ 


(3.8) 


One of the reasons of the interest of our results, in particular of the first-integrals (13.61) and (13.8L is that 
they provide a new geometric approach to results provided previously in [^. Additionally, along with 
Proposition [521 it enables us to obtain new Lagrangian descriptions of the oscillators under study. 

Let us work out a Lagrangian for (13.2p via /ii. From equation 


d'^Li 

dv'^ 


= = \ + kx^, 


we obtain 


Li{x,v) = -(1 -I- kx^)v‘^ + (l>i{x)v + 4’2 {x), 


for certain functions (j)i and 4>2- The gauge term (j>i{x)v can be set to zero and 02(a^) is determined by 
p.lll) where we choose uq = 0, i.e. 


02(a;) = [ (/xF)(C,0)dC = [ 

Jo Jo 


0 (l + fcC^)^ 


dC = 


a 


2k{l + kx^) 


2 2 
a 


2(1 -I- kx^) ' 


Then, the Lagrangian for (13.21) is given, up to addition of a gauge term, by 


Li{x,v) 


i(l -I- kx^) 


1 x'^ 

2 1 + kx'^ ’ 


with corresponding Hamiltonian 

Hi{x,p) 


1 1 x^ 

2 1 - 1 - kx"^ ”*"2 1 - 1 - kx'^ 


Meanwhile, if we make use of the second Jacobi multiplier, p ,2 = ^(1 + kx’^Y'v^ — we obtain a 
non-standard Lagrangian, because 


d‘^L2 

dv'^ 


= P2^ L2{x, v) = + kx^f - 


02 (a;). 
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with 


M-} = 1 (Mm.o)d< = ^ 

Hence, up to an irrelevant gauge term, we obtain that 


ax‘^{2 + kx^) 
4(1 + kx'^Y 


L2{x,v) 


^ 4/1 I 7 2\2 ^ ^ 

-V (1 + kx) - — 


ax ^(2 + kx'^) 
4(1 + kx'^)^ 


Consequently, 


k , 2sO -i 2 , \/p+^ 

p = —(1 + kx^Yv^ — a^v => u = ± , ^ —7 

3 Yk + 2k'^x"^ + k^x'^ 


and 


H2{x,p) = 


p^ + 4a^ + a{—kx‘^{2 + kx"^) + 3a^ — Aa^Jk{\ + kx‘^Y{p + a^) 

4fc(l + /cx^)^ 


4 Jet bundle formulation of symmetries of differential equations 

Systems of differential equations and their symmetries admit an alternative geometric approach: instead 
of considering them as vector fields, we describe them as geometric structures in jet bundles. We now 
recall the basic ingredients of this formulation before passing to study non-local symmetries of differential 
equations. 

Let {E, K., tt) be a fiber bundle with total space £1 = R x iV, base R and submersion tt : (f, x) G 
E t G M.. Given local coordinates {x^}j=i,...,n on N and t on M, we can naturally define a coor¬ 
dinate system = x^}j=i^,,,^n on E. Given a section cr : R —>■ i? of (£1,R, 7 r), we write = 

{to,xopXip... ,Xk)) for the fc-order jet prolongation of cr at to, he. the equivalence class of sections 
7 : t S R !->• (t, 7 ^(t),... , 7 ”(t)) G E such that 

Y{to)=x^Q), -^{to)=xl-^, i = l,...,fe, j = l,...,n. 

We denote by J^tt, for £ > 0, the space of £-order jet prolongations (fc-jets) of the fiber bundle (£',R, tt) 
and we define J° 7 r = E. Alternatively, we write J^(R, A) for J^’tt when tt is understood from the knowl¬ 
edge of E and R. The space is a finite-dimensional manifold with local coordinates {t,x^.,} i=o. k 

of the form 


i = 0,...,k, j = l,...,n, 

with d°7'l/dt°(to) = 7'^(to)- It is well known that (J^7 r,R,: J^tt —>■ R) is a fiber bundle: the k-order 
jet bundle associated with (£1, R, tt ). The sections of the fc-order jet bundle being the prolongation of a 
section of E are called (fc-order) holonomic sections l36\ 137] . 

Gonsider the C°°(J^ 7 r)-module of 1-forms 6 G A^(J^ 7 r) satisfying that {j^<j)*0 = 0 for every section 
tr : R —>■ £1. The elements of this module are called contact forms on J^tt. It can be proved that this 
module is a locally free-module generated by the contact forms 0 ^^ = dx^^ — x^_|_j^j dt, with j = 1 ,... ,n 
and i = 0,..., fc — 1. We can endow J^tt with the distribution spanned by vector fields annihilating 
all contact forms on J^tt, the referred to as contact or Cartan distribution of J^tt. In particular, tangent 
vectors to graphs of fc-order jet prolongations belong to C^. More generally, the Gartan distribution of 
J^TT is spanned by 


n k—1 

= ^ + X] X] ^i+l) 

j=l i=0 




dxj^ 


D, = 


dxi 


j = 1 , 


k) 


(4.1) 
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Note that is not involutive for finite k and have dimension n + 1. 

We call Lie symmetries'^ of the infinitesimal symmetries of i.e. the vector fields Y on J^-k 
satisfying that takes values in for every vector field X taking values in C^. In other words, 

Lie symmetries of C* are those vector fields whose flows give rise to transformations mapping /c-order 
holonomic sections into fe-order holonomic sections. 

Given a vector field X on E, its prolongation to J^tt is the unique Lie symmetry, X^^\ of whose 
holonomic integral curves are the prolongations to J^tt of integral curves of X. Equivalently, X^^'> is the 
unique Lie symmetry of projecting onto X under Tikp '■ jt<^ G J^tt i—> cr(t) G E. Lie point symmetries 
of are Lie symmetries of that are prolongations to of vector fields on E. The Lie-Backlund 
theorem states that not all Lie symmetries of are Lie point symmetries. Non-Lie point symmetries, 
the referred to as contact Lie symmetries, can always be considered as liftings of a uniquely defined Lie 
symmetry on J^tt (see m)- 

In the above framework, a fc-order system of differential equations is a closed embedded submanifold 
£ C and its particular solutions are sections of tt whose /c-order prolongations belong to £. We say 
that a system of k-order differential equations is in normal form when the natural projection iTk^k-i ■ 
G £ ^ G J^~^TT is a submersion, e.g. £ = {jfx G J^(R,K^) : ~ * 2 ) = 

say that £ is in normal form and not underdetermined when iTk^k-i ■ G £ G is 

a diffeomorphism, e.g. £ = {j^x G J^(R,R^) | x^^ = E^{t,x,xi)),j = 1,2} for arbitrary functions 
F^,E'^ : J^(R,R^) —>■ R. When £ C J^tt is in normal form and not underdetermined, there exists a 
vectorial mapping A : J^tt R" allowing us to write that £ = A“^(0). 

A Lie symmetry of that is tangent to £ is called a classical symmetry of £ |15j . On the other hand, 
the term classical symmetry of £ has also being employed |14| to refer to a vector field on E giving rise 
to a uniparametric group of transformations mapping particular solutions to £ to particular solutions to 
£. In this work we will mainly use the first definition. Nevertheless, if a classical symmetry E is a Lie 
point symmetry, then we also call classical symmetry the unique vector field on E whose prolongation to 
J^TT is Y. More specifically, we say that E is a classical point symmetry of £. 

The projections TTh,k '■ J^tt, with h > k, enable us to define the bundle of infinite jets 

(J°°7r, R, TToo : J°°TT —?► R) as the inverse limit of the projections 

R ^ E ^ J^TT ^ J^TT ^ . 

The commutative ring of differentiable functions over J°°tt is defined by E{tt) = [j’^QC°°{j’‘TT). Note 
that each element of depends on a finite subset of variables of {t,xL} igNu{o} ■ When there exists a 

j=l,....n 

natural injection between two manifolds, e.g. J^tt ^ tt for k' > k, we may consider each function on 
the first manifold as a function on the second, e.g. an element of C°°{J^'it) as an element of tt), 

so as to simplify the notation. 

Given a section cr : R —>■ J^tt, we call infinite prolongation of a the section j°°a : t G R 1 —>• jf°cr G J°°tt 
given in coordinates by 

j“cr= . 

Similarly to finite-dimensional manifolds, vector fields on J°°tt are defined as derivations of the com¬ 
mutative ring J-{tt) and the J'(7r)-module of vector fields on J°°tt becomes a Lie algebra with respect 
to the commutator of derivations. The tangent vectors to infinite prolongations of sections of E span a 
distribution on J°°tt spanned by the derivation on E{tt) given by 

Tl OG Q 

Although D depends on an infinite number of variables, it induces a well-defined derivation on E{tt) due 
to the fact that every function in E{tt) only depends on a finite number of variables. The distribution C 
spanned by D is the referred to as contact or Cartan distribution of J°°tt, which is one-dimensional and 
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therefore involutive. We cannot ensure that D is integrable as the Frobenius Theorem does not apply to 
distributions on infinite-dimensional manifolds. 

We call Lie symmetries of C the infinitesimal symmetries of C. As these Lie symmetries of C are 
defined on an infinite-dimensional manifold, we cannot ensure that they are associated to uni-parametric 
groups of diffeomorphisms on 

It is interesting to note that many structures on the jet spaces J^tt become simpler when passing to 
e.g. the Cartan distribution becomes one-dimensional and involutive. Moreover, J°°7r is geometri¬ 
cally richer than finite-order jet bundles. For instance, Lie symmetries of C need not be lifts neither of 
vector helds on E nor of vector fields on J^tt due to the fact that the Lic-Backlund theorem does not 
apply to such Lie symmetries m- 

Given a fc-order system of differential equations £ C J^tt, the I-prolongation of £ is the set of points 
£(1) = {jf+‘a\fa C£}C (see for details). Further, we call infinite prolongation of £ the 

set £°° = I j^cr C £} C J°°Tr. If £ is in normal form and not underdetermined, then £ = A“^(0) 
for a certain mapping A : jfx G J^tt >->• ... ,x^^ — G K” and functions 

Fj ■ case, £°° is a finite-dimensional manifold locally determined by the infinite set 

of conditions 

s—times 

A = 0, D^A=D{D{...D{D A)...)) = 0, (4.2) 

with s = 1,2,... The above conditions determine all derivatives of particular solutions of £ out of the 
value of the first {k— l)-derivatives. Hence, a local set of coordinates for £ can be considered in a natural 
way as elements of Fljr) giving rise to a coordinate system on £°^, which becomes a finite-dimensional 
manifold. In view of (14.2L the restriction of D to £°° is tangent to This allows us to endow 
with a one-dimensional distribution C\s°°- The pair given by £^ and C\s°° becomes what is called a 
one-dimensional diffiety. We define F{£°°) to be the restriction to £°° of functions of .F(7r). The Lie 
symmetries of C that are tangent to £°° are called higher symmetries of £. 

Definition 4.1. Let f be a fc-order system of ODEs, we say that a bundle {£^,£°°,k : £'^ —>■ £°°) is a 
covering for £ if the bundle £‘^ can be endowed with a one-dimensional distribution 

= {Cp}pe? 

in such a way that (k*)^ : Cf —?► C«;(p) is a linear isomorphism for each p G £‘^. 

Definition 4.2. Let f be a fc-order system of ODEs and let £'^ be a covering for We call non-local 
symmetry for £ an infinitesimal symmetry of the distribution C^. 

We call dimension of the covering, dim(K), the dimension of the fiber of the bundle k. Given a covering 
K, integral manifolds S of C‘^ project under k onto integral manifolds of Cjgoo, i.e. onto prolongations of 
particular solutions of £. It follows that, in this picture, symmetries of shuffle integral manifolds of 
and, projecting under k, we can obtain particular solutions to £. 

5 Extended formalism and non-local symmetries 

Recently Gandarias and coworkers have studied a method for obtaining non-local symmetries for certain 
particular systems of ODEs [SJEllT]. The main idea is to add a new additional equation to the original 
one in such a way that the new higher-dimensional system can be endowed with a classical symmetry. 
In this section, we study this technique from a geometric perspective. The starting point is that this 
procedure must be considered, in geometric terms, as a very particular case of an extended formalism 
where the covering has a total space given by a one-dimensional diffiety. We also mention the relation of 
this approach with some results previously obtained by Krasil’shchi, Vinogradov and coworkers dsiiisiis]. 

We begin with a system £ of fc-order ODEs given by £ = A~^(0) C for a mapping A : 

J^Tn+i —>■ RP with (R"+^, R, r„+i : (t,x) G R x R" i—>■ t G R) and coordinates 

/ H T H ^ T \ 

AM<,x,—=0, i = l,...,p, xgR" (5.1) 
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and we assume that E is underdetermined and in normal form, i.e. n = p and iTk^k-i ■ £ C J^tt —>■ 
J^~^TT is a diffeomorphism. A classical point symmetry of f, represented by an e-parametric group of 
transformations given infinitesimally by 

1-7 7 7 \ ^ 1,..., n, 

lx-' = x-' + erp [t, X), 


preserves the set of solutions of the equation, that is, transforms particular solutions of (EH) into particular 
solutions of the same equation. Unfortunately, many systems of differential equations do not possess 
classical point symmetries. In that case, the method of non-local symmetries (applied to some particular 
cases in ISEIIT] and related to some questions discussed in [El [mis]) can be of a great usefulness. 

A classical point symmetry for the system EH is a vector field Y on such that (i) is the lift 

of a vector field Yq on ~ and (ii) is tangent to the submanifold £. In coordinates if Yq 

takes the form 

Yoit.x) = —, 


then Y is given by 




dt 


+E 

i=i 






(5.2) 


where the functions ^ : J Tn+i C J°°Tn+i —t R can be obtained from the functions ^ and , with 
j = 1,..., n (for a detailed explanation see|38]l. 

Suppose that the system EH is given. Let us construct a new system containing EH as a particular 
part. Consider the new jet bundles associated to the bundle (M”+^,R, t „+2 : (t,x) € R”+^ ^ t G R), 
with X = {x,w) G R"+^. For a certain fixed function H : J^Tn+i —t R, we can define a new system 
£ C J^Tn +2 as follows 


. / dx 




w 


dt2 


= Dk)H. 


dt^ 


= 


(5.3) 


containing as a particular part the initial system EH and where Dk) is one of the generators of the Cartan 
distribution of J^r „+2 given in (14.11) . Observe that although £ was in normal form, the new system (15.31) is 
not in normal form as well: the derivatives d^w/dt^ can be expressed in terms of the lower derivatives, but 
such lower derivatives must hold several additional conditions. More geometrically, we cannot consider 
(15.31) straightforwardly as a submanifold £ C J^Tn +2 in such a way that Tk._k-i ■ £ C J^Tn +2 —t J^~^Tn +2 
is epijective. ^ ^ 

In any case, we can consider the prolongations £°^ and of £ and £, respectively. Since the system 
£ is in normal form and in view of the dehnition of £, we see that the values of each derivative d^x/dt^, 
with p > k, and d^w/dtP, with p > 1, of each particular solution of £ and £ can be determined from the 
previous derivatives. Thus, and £°° are finite-dimensional manifolds and 

dim f = dim £), dim£)°“ = dimf. 


Hence, a local coordinate system on £ or £ induces a local TOordinate system on their infinite prolonga¬ 
tions. This means that expressions in coordinates on £ and £ can be understood as expressions on £°° or 
£°° indistinctly. This property is important: it allows us to identify £°° with £ and £^ with £. Hence, 
calculations in infinite-dimensional jet bundles are as difficult as for hnite-dimensional jet bundles and 
the whole procedure is properly defined in a rigorous more powerful geometrical manner. For instance, 
D, which has no sense on j’^Tn+i, can be however correctly considered when restricted to £°°. 

Since D is tangent to we can define the restriction D\£^ of this operator to If D is the 
analogue of D on J°°t„+ 2 , this operator is also tangent to £°° and we can also define its restriction, 
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to S°°. The vector field D induces a one-dimensional distribution C on J°“t „+2 and spans a 
one-dimensional distribution turning the pair {£°°into what is called a diffiety of dimension 

one: the dimension one refers to the fact that we have defined a one-dimensional distribution on 

£°°. Moreover, we have the following property. 

Proposition 5.1. Let D he the vector field on J°°Tn +2 given by 


D = D 


d 


OO 

i =0 




Wo = w. 


(5.4) 


Then, we have locally 


T’Icoo = D\eo 


a n 

‘'l-T. 


‘k-2 

T.4+1) 

.p—O 


H) 


+ F^Ut^) 


M-i) 




(5.5) 


for certain functions F^,..., F" : £°° —>■ R. 

Proof. Expression (EH) is trivially the generator of the Cartan distribution of J°°t„+ 2 . As we assume £ 
to be in normal form and not underdetermined, the higher-order derivatives of the particular solutions 
to £, namely d^a;'^7dt^, can be locally written as a function of the previous derivatives. So, d^a;'^7dt^’ = 
F ^for j = 1,..., n and certain functions F^ : —>■ K which is understood in the natural 

way as a function on £^ C J°°Tn+ 2 - Using this, recalling that {t,Xi),w}i=o,...,k-i forms a coordinate 
system for £°° and restricting D from J°°Tn +2 to £°°, we obtain that the expression (15.51) follows from 
( 1 ^ . □ 

The natural projection If : {t, x, w) € R"+^ i->- (t, x) G lifts to a projection J°°n : J°°t „+2 —>• 

J°°Tn+i satisfying 

j^uurx) = j°^n{jr{x,w)) =jrx. 

This projection induces a map : £°^ —>• £°° obeying that 


As a consequence, J°°n induces an isomorphism (J°°n7oo)*j : Q — >• Cjoon( 5 ) for every f G £°° and 
therefore a covering k* = (J°“n|goo), : C\s°°. Hence, if the system £ admits a classical symmetry 

Y, e.g. (15.2L then Y can be lift to a Lie symmetry Y°^ of C, namely a higher symmetry for £. Both vector 
fields, Y and Y°°, are tangent to £ and respectively. It is worth noting that due to our assumptions, 
the coordinate expression in £ of Y\^ and the coordinate expression of in £°° are the same. 

Moreover, as is a higher symmetry, it leaves invariant C and it becomes a non-local symmetry of £'^ 
when restricted to £°°. 

If £ admits a classical point symmetry Y, then we have a one-parametric group of diffeomorphisms 
given, in infinitesimal form, as 

{ t* =t + ef{t,x,w ), 

X* = a; -I- e (j){t, x, w ), 
w* = w + e rj{t, X, w ), 

transforming solutions to the system £ into solutions of £. Hence, the set of transformations 


t* =t + ef{t,x,w{f )), 
a;* = a; -|- e 4>(t, x, w(t )), 


enables us to map particular solutions to £ into solutions of £ by means of the curves w{t) corresponding 
to particular solutions of £. 
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We can summarise the main results proved in this section as follows. We have proved we can embed 
a given system E into a bigger one whose infinity prolongation has the structure of a difSety. This 
structure gives rise to a covering for the initial system. In |15| another covering for the initial system 
is constructed so as to study it through non-local symmetries. Nevertheless, that covering cannot be 
considered neither as a diffiety nor a submanifold of a jet bundle without additional consteuctions. So, 
our interpretation is more powerful as we can straightforwardly use classical symmetries of £ to construct 
non-local symmetries of £. Gandarias developed slightly modifications of her method, but all of them 
can be retrieved as particular cases of the above geometric approach. 

6 Extended formalism for the nonlinear oscillators 

In this section, we provide a geometric method to construct non-local symmetries for second-order au¬ 
tonomous differential equations. This method is based upon considering our initial second-order differ¬ 
ential equations as part of an extended system whose form can be determined out of the initial one. 
In following subsections we show that this procedure retrieves as a particular case the results given by 
Gandarias and coworkers. 

Theorem 6.1. Every system £ corresponding to 

^^ = F{x,v), (6.1) 

can he extended to a larger system £ with the additional equation dw/dt = H{x,v) in such a way that 
Y = gXfj, where Xjj = vd^ + Edy + Hd^j is an infinitesimal symmetry of the distribution generated 
by dt + Xh and g G J-(£°°) is a first-integral of Xh, dx, dt (as vector fields on £°°). In consequence, Y 
is a non-local symmetry of (EJP. 

Proof. We have that [gXH,dt d- Xh] = —{dtg + XHg)XH. So, gXn is an infinitesimal symmetry of 
C|goo if and only if dtg Xng = 0. Under the assumed conditions for g, we obtain that gXn is an 
infinitesimal symmetry of Let us prove that there exists a nonconstant function g satisfying the 

above conditions. 

Set H{x,v) = F{x,v)h{v) for a certain function h{v) 0. Since dgjdx = 0, then 

li F 0 and we require g to be non-constant, then our definition of H ensures the existence of a non¬ 
trivial g depending only on v and w. li F = 0, we can choose any g with the required properties of our 
theorem. □ 

Corollary 6.2. A classical infinitesimal symmetry for the system £ gives rise to a non-local symmetry 
of£. 

Proof. Every classi£al symmetry for £ can be extended in view of the Lie-Backlund Theorem to a higher 
symmetry Y°^ of £. This higher symmetry is tangent to £°° and a Lie symmetry of C, which is, by 
construction of tangent to £°°. Hence, Y°°\g„^ is a symmetry of C\g^ and it becomes a non-local 
symmetry for £. □ 
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6.1 First nonlinear oscillator 


Let us review the approach given by Gandarias to study the nonlinear oscillator (I3.1|) . The first-order 
system (13.31) associated to (13.11) is embedded into a new one on TR x K of the form 


dx 

dn — a^)a 


dt 1-1- kx"^ 


( 6 . 2 ) 


where H(x,v) is a, undetermined for the moment, function, to be fixed later on. Let us study the Lie 
point symmetries of this system. Particular solutions to system (16.21) are in a one-to-one correspondence 
with the integral curves t >->• {t, x{t),v{t),w{t)) of the vector field 




dx 


— a‘^)x d TTf \ ^ 

l + kx^ ^ 


Given a vector field Y = ^d/dt + (l>d/dx + iljd/dv + r]d/dw on x TR, where we include the time variable 
t, we know that Y determines a Lie point symmetry of this system if [T, dt + Xh] = f{dt + Xh), where 
we recall that t,x,v,w are considered as coordinates on R^ x TR and / G x TR). Equivalently 

E is a classical point symmetry of this system if its four coefficients satisfy the following equations 


A^)a=o = v(kx^ + 1)6 + vH{kx^ + 1).^^ -|- xv{kv‘^ - a^)6 + v^{kx^ + 1 )^ 0 , 

— x{kv^ — a‘^)4>v — ikx^ + l)v(j)x — {kx^ + l)(f)t — H{kx^ + l)(j)w + v^{kx^ + l)^/i = 0, 

— Hxk(j)x^ -I- H{kx^ -I- l)r]yj + rjtkx^ — a^xrjy — Hytp — Hx(l> -f ?7x = 0, 
(y(i)A^)A=o = {kx"^ + l)(a!^ — kv‘^)x^t — H{kx^ + \)x[kv^ — a^]^w — x‘^{kv'^ — a^)^6 

-I- vx[l + kx^]{—kv‘^ + a^)^x + {kx^ + l){kv^ — a^)xil}y + {kx^ — \){kv‘^ — a^)4> 

+ {kx^ + l)'^vipx — 2kxv{l + kx^)il} + H(kx^ + l)^'0iu + {kx^ + l)^<)>t = 0 
(y(i)A3)A=o = -H{kx^ + 1)6 - H^{kx^ + l)e^ - Hx{kv^ - a^)^y 

— H[kx^ + l)v^x — Hx^[kx^ + 1 ) + kxv'^r]y + {kx^ + l)vr]x — Hyk^jjx^ 

(6.3) 


for A^ = i; — u, A^ = u — {kv^ — oP’^xj (1 -I- kx'^), = w — H and Y^'^ being the prolongation to J^r 4 , 
with T 4 : {t,x,v,w) € R^ >->• t € R, of the vector field Y on R^. We include expressions (16.3() to solve 
several minor typos and mistakes in the previous literature. This is a quite difficult system to be solved, 
which suggests us to assume some kind of simplification. This was done in [5] , whose authors considered 
as an ansatz a particular form for 6 <(', I- Now we reconsider this whole approach in a more geometrical 
and rigorous way. 

Equivalently, the differential equation dUD can be considered along with the equation dw/dt = 
H{x,v). As commented in the latter section, this system can be understood as a submanifold £ of 
with ra : {t,x,w) £ R^ i—>■ t G R. Let us use Theorem 16.II to study the infinitesimal symmetries of Cjgoo. 
Recall that this amounts to a non-local symmetry for £. 

We can construct a non-local symmetry by assuming Y = gXn with H{x,v) = {XHv){x,v)h{v) and 
g being a first-integral of Xh independent of t and x, namely, such that 


Xng 


{kv'^ — a^)x 
1 -I- kx'^ 


dg ^ dg 


= 0 , 


where we fixed according to Theorem 16.II 


H{x, v) 


{kv^ — a^)x 
1 -I- kx'^ 


h{v). 
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By assuming h{v) = 1/v, we obtain a simple first-integral for Xh of the form g = jv. Hence, 


F = e“ 


A nr/ A 

dx dv 


H d 
V dw 


Indeed, observe that \Y,Xh\ = 0. 

As {t^x,v,w} can be understood as coordinates of x TR and S, the vector field Y can also be 
considered as a vector field on R^ x TR. In this way, Y is the same Lie symmetry provided in [^, where 
it was obtained by the derivation of a particular solution of (16.31) using an ad hoc ansatz for Y and H. 
Meanwhile, we here use a covering to show that Gandarias’ and coworkers ansatz corresponds to choose 
a certain H so that a Hrst-integral for Xh independent of a;, t can be obtained. This immediately leads 
to their same hnal result. 

Note also that we could in principle choose another function H which could potentially lead to different 
non-local symmetries of £. Nevertheless, the form chosen in this work makes computations easier in many 
cases. 


6.2 Second nonlinear oscillator 

We can now apply the above method to equations (13.21) to recover the same result provided in . In this 
new case, the first-order system (Irfl) associated to (13.21) is embedded into one 


da: 

dt 

du kxv^ 
dt 1-1- kx'^ 
dw 


(1 -I- kx'^y ’ 


dt 


= H{x,v), 


(6.4) 


on R^, where H{x,v) is a function to be fixed next. Additionally, we can consider this system as a 
submanifold £ C with ra : {t,x,w) G R^ >->• t € R. This system describes the integral curves 
t —(t, a:(t), v{t),w{t)) of the vector field on R^ x TR ~ f of the form 


- _ 9 ^ ^ f kxv'^ o?x \ 9 9 

^ 9t ^ ^ dx \1 + kx'^ {1 + kx'^Y J dv dw 

We fix H to be of the previously commented form, i.e. 

/ kxv"^ a^x \ 1 

{X,V) = - + j -■ 


Hence, Xh admits a locally dehned first-integral g that does not depend neither on x nor on 
such that 




/ kxv^ o?x \ dg dg 

\1 -I- kx^ (1 -I- kx'^Y ) dv dw 


namely, 


This leads to a simple first-integral for Xh of the form g = jv. We can now obtain a Lie symmetry of 
the system by choosing Y = gXn, which reads 


dx dv V dw 


which is again the same classical symmetry provided in but we here understand it as a symmetry of 
C on £°^, i.e. a non-local symmetry of £. 


17 











7 Final comments 


This paper has been mainly concerned with the following two points: Jacobi multipliers and non-local 
symmetries. 

• The Jacobi multipliers have been first considered in relation with the inverse problem of the La- 
grangian formalism and then applied to the study of two particular nonlinear oscillators. 

• The theory of non-local symmetries has been studied by making use of a geometric approach. We 
prove that the extended formalism can be a very interesting procedure for obtaining symmetries of 
nonlinear systems. 

• We have shown that the use of infinite-dimensional jet manifolds does not complicate the description 
of non-local symmetries of systems and allows us to develop a more rigorous theoretical approach. 
In addition, certain structures are now naturally defined. 

• In the future we aim to apply the theory of non-local symmetries to a generalisation of the nonlinear 
oscillators studied in this work that contain an isotopic term. This will describe as a particular case 
the non-linear oscillators detailed in [10] on a one-dimensional manifold. 

• Diffieties are mainly used in the study of systems of partial differential equations. Nevertheless, we 
aim to show that these structures may play a role also for the study of relevant systems of first-order 
differential equations. 
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